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Modified gravity theories generically predict a violation of Lorentz invariance, which may lead 
to a modified dispersion relation for propagating modes of gravitational waves. We construct a 
parametrized dispersion relation that can reproduce a range of known Lorentz-violating predictions 
and investigate their impact on the propagation of gravitational waves. A modified dispersion 
relation forces different wavelengths of the gravitational wave train to travel at slightly different 
velocities, leading to a modified phase evolution observed at a gravitational-wave detector. We 
show how such corrections map to the waveform observable and to the parametrized post-Einstcinian 
framework, proposed to model a range of deviations from General Relativity. Given a gravitational- 
wave detection, the lack of evidence for such corrections could then be used to place a constraint 
on Lorentz violation. The constraints we obtain are tightest for dispersion relations that scale with 
small power of the graviton's momentum and deteriorate for a steeper scaling. 
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I. INTRODUCTION 

After a century of experimental success, Einstein's fun- 
damental theories, ie. the special theory of relativity and 
the General theory of Relativity (GR), are beginning to 
be questioned. As an example, consider the observation 
of ultra-high-energy cosmic rays. In relativity, there is a 
threshold of ~ 5 x 10 19 eV (GZK limit) for the amount of 
energy that charged particles can carry, while cosmic rays 
have been detected with higher energies pQ . On the theo- 
retical front, theories of quantum gravity also generically 
predict a deviation from Einstein's theory at sufficiently 
large energies or small scales. In particular, Lorentz vi- 
olation seems ubiquitous in such theories. These con- 
siderations motivate us to study the effects of Lorentz 
violation on gravitational wave observables. 

Einstein's theory will soon be put to the test through 
a new type of observation: gravitational waves (GWs). 
Such waves are (far-field) oscillations of spacetime that 
encode invaluable and detailed information about the 
source that produced them. For example, the inspiral, 
merger and ringdown of compact objects (black holes or 
neutron stars) are expected to produce detectable waves 
that will access horizon-scale curvatures and energies. 
Gravitational waves may thus provide new hints as to 
whether Einstein's theory remains valid in this previously 
untested regime. 

Gravitational-wave detectors are today a real- 
ity. Ground-based interferometers, such as the Ad- 
vanced Laser Interferometer Gravitational Observatory 
(Ad. LIGO) [5H3] and Advanced Virgo [5], are currently 
being updated, and are scheduled to begin data acquisi- 
tion by 2015. Second generation detectors, such as the 
Einstein Telescope (ET) [U [7] and the Laser Interferom- 
eter Space Antenna (LISA) 8 , 9 J , are also being planned 



for the next decade. Recent budgetary constraints in the 
United States have cast doubt on the status of LISA, 
but the European Space Agency is still considering a 
descoped, LISA-like mission (an NGO, or New Grav- 
itational Observatory). The detection of gravitational 
waves is, of course, not a certainty, as the astrophysical 
event rate is highly uncertain. However, there is consen- 
sus that advanced ground detectors should observe a few 
gravitational- wave events by the end of this decade. 

Some alternative gravity theories endow the graviton 
with a mass [10] . Massive gravitons would travel slower 
than the speed of light, but most importantly, their speed 
would depend on their energy or wavelength. Since grav- 
itational waves emitted by compact binary inspirals chirp 
in frequency, gravitons emitted in the early inspiral will 
travel more slowly than those emitted close to merger, 
leading to a frequency-dependent gravitational-wave de- 
phasing, compared to the phasing of a massless gen- 
eral relativistic graviton. If such a dephasing is not ob- 
served, then one could place a constraint on the graviton 
mass |11) . A Lorentz-violating graviton dispersion rela- 
tion leaves an additional imprint on the propagation of 
gravitational waves, irrespective of the generation mech- 
anism. Thus a bound on the dephasing effect could also 
bound the degree of Lorentz violation. 

In this paper, we construct a framework to study the 
impact of a Lorentz-violating dispersion relation on the 
propagation of gravitational waves. We begin by propos- 
ing a generic, but quantum-gravitational inspired, mod- 
ified dispersion relation, given by 

E 2 =p 2 <? +m 2 c 4 + Ap a c a , (1) 

where m g is the mass of the graviton and A and a are 
two Lorentz-violating parameters that characterize the 
GR deviation (a is dimensionless while A has dimensions 
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of [energy] 2- "). We will assume that A/(cp) 2 ~ Q <C 1. 
When either A = or a = 0, the modification reduces to 
that of a massive graviton. When a = (3, 4), one recovers 
predictions of certain quantum-gravitation inspired mod- 
els. This modified dispersion relation introduces Lorentz- 
violating deviations in a continuous way, such that when 
the parameter A is taken to zero, the dispersion relation 
reduces to that of a simple massive graviton. 

The dispersion relation of Eq. modifies the gravi- 
tational waveform observed at a detector by correcting 
the phase with certain frequency-dependent terms. In 
the stationary-phase approximation (SPA), the Fourier 
transform of the waveform is corrected by a term of 
the form £(A)w Q-1 , where u — irAif is a dimension- 
less measure of the gravitational-wave frequency with M 
the so-called "chirp mass". We show that such a mod- 
ification can be easily mapped to the recently proposed 
parametrized post-Einsteinian framework (ppE) [121 113] 
for an appropriate choice of ppE parameters. 

In deriving the gravitational-wave Fourier transform 
we must assume a functional form for the waveform as 
emitted at the source so as to relate the time of arrival at 
the detector to the gravitational-wave frequency. In prin- 
ciple, this would require a prediction for the equations of 
motion and gravitational-wave emission for each Lorentz 
violating theory under study. However few such theories 
have reached a sufficient state of development to produce 
such predictions. On the other hand, it is reasonable to 
assume that the predictions will be not too different from 
those of general relativity. For example, we argued 
that for a theory with a massive graviton, the differences 
would be of order (A/A g ) 2 , where A is the gravitational 
wavelength, and X g is the graviton Compton wavelength, 
and A 9 > A for sources of interest. Similar behavior 
might be expected in Lorentz violating theories. The im- 
portant phenomenon is the accumulation of dephasing 
over the enormous propagation distances from source to 
detector, not the small differences in the source behavior. 
As a result, we will use the standard general relativistic 
wave generation framework for the source waveform. 

With this new waveform model, we then carry out a 
simplified (angle-averaged) Fisher-matrix analysis to es- 
timate the accuracy to which the parameter C(-^0 could 
be constrained as a function of a, given a gravitational- 
wave detection consistent with general relativity. We per- 
form this study with a waveform model that represents 
a non-spinning, quasi-circular, compact binary inspiral, 
but that deviates from general relativity only through the 
effect of the modified dispersion relation on the propaga- 
tion speed of the waves, via Eq. 0. 

To illustrate our results, we show in Table [I] the accu- 
racy to which Lorentz-violation in the a = 3 case could 
be constrained, as a function of system masses and detec- 
tors for fixed signal-to- noise ratio (SNR). The case a = 3 
is a prediction of "doubly special relativity" . The bounds 
on the graviton mass are consistent with previous stud- 
ies [111 I14H18) (for a recent summary of current and pro- 
posed bounds on m g see [IH])- The table here means that 



given a gravitational-wave detection consistent with GR, 
m g and A would have to be smaller than the numbers on 
the third and fourth columns respectively. 



Detector 
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A(eV" 1 ) 


Ad. LIGO 


1.4 


1.4 


3.71 x 10" 




7.36 x 10" 8 


SNR = 10 


1.4 


10 


3.56 x 10" 




3.54 x 10" 7 




1 n 


1U 


^ m v i 

O.dl A 1U 


22 


D.oo X 1U 


ET 


10 


10 


2.99 x 10" 


23 


2.32 x 10" 8 


SNR = 50 


10 


100 


4.81 x 10" 


23 


1.12 x 10" 6 




100 


100 


6.67 x 10" 


23 


3.34 x 10" 6 


NGO 


10 4 


10 4 


3.05 x 10" 


25 


2.16 x 10" 2 


SNR = 100 


1() 4 


10 5 


2.46 x 10" 


25 


0.147 




10 s 


10 5 


2.03 x 10" 


25 


0.189 




10 s 


10 6 


2.09 x 10" 


25 


9.57 




10 6 


10 6 


1.49 x 10" 


25 


23.2 



TABLE I: Accuracy to which graviton mass and the Lorentz- 
violating parameter A could be constrained for the a = 3 case, 
given a gravitational-wave detection consistent with GR. The 
first column lists the masses of the objects considered, the 
instrument analyzed and the signal-to-noise ratio (SNR). 



The remainder of this paper deals with the details of 
the calculations and is organized as follows. In Sec. [TTJ 
we introduce and motivate the modified dispersion rela- 
tion ([I]), and derive from it the gravitational- wave speed 
as a function of energy and the new Lorentz-violating 
parameters. In Sec. |III[ we study the propagation of 
gravitons in a cosmological background as determined by 
the modified dispersion relation and graviton speed. We 
find the relation between emission and arrival times of 
the gravitational waves, which then allows us in Sec |IV| 
to construct a restricted post-Newtonian (PN) gravita- 
tional waveform to 3.5 PN order in the phase [0(v/c) 7 }. 
We also discuss the connection to the ppE framework. 
In SecjVj we calculate the Fisher information matrix for 
Ad. LIGO, ET and a LISA-like mission and determine 
the accuracy to which the compact binary's parameters 
can be measured, including a bound on the graviton and 
Lorentz-violating Compton wavelengths. In Sec. |VI| we 
present some conclusions and discuss possible avenues for 
future research. 



II. THE SPEED OF LORENTZ-VIOLATING 
GRAVITATIONAL WAVES 

In general relativity, gravitational waves travel at the 
speed of light c because the gauge boson associated with 
gravity, the graviton, is massless. Modified gravity theo- 
ries, however, predict modifications to the gravitational- 
wave dispersion relation, which would in turn force the 
waves to travel at speeds different than c. The most 
intuitive, yet purely phenomenological modification one 
might expect is to introduce a mass for the graviton, fol- 
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lowing the special relativistic relation 

1^2 2 2 i 2 4 

E = p c + mtc ■ 



(2) 



From this dispersion relation, together with the defini- 



tion v/c = p/p°, 
speed [IT] 



or v — c 2 p/E , one finds the graviton 



= 1 - 



2 4 
m 9 C 

E 2 ' 



(3) 



where m g , v g and E are the graviton's rest mass, velocity 
and energy. 

Different alternative gravity theories may predict dif- 
ferent dispersion relations from Eq. A few examples 
of such relations include the following: 

• Double Special Relativity Theory 20-23 : E 2 — 
p 2 c 2 + m 2 c 4 + i] dslt E 3 + where ?y dsrt is a pa- 
rameter of the order of the Planck length. 
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• Extra- Dimensional Theories [24j : E 2 
m 2 c 4 — ttedt-E 4 , where a cdt is a constant related 
to the square of the Planck length; 



• Hofava-Lifshitz Theory [25H27] : E 2 = p 2 c 2 + 

( K hiMhiA6) P 4 + • ■ •' where K hi an d Hm are con- 
stants of the theory; 

• Theories with Non-Commutative Geometries |28j : 
E 2 gl{E) = m 2 g c i +p 2 c 2 gl{E) with g 2 = 1 and 
.91 = (1 - yta ncg Tr/2)eTip(-a ncg E 2 /E 2 ), with a ncg 
a constant. 



rrtg , but small compared to the Planck energy E p . More 
precisely, we have 

• Double Special Relativity: A = ?7d S rt and a = 3. 

• Extra-Dim. Theories: A = — a c A t and a = 4. 



Hofava-Lifshitz: 



with 



Mil 



£*hi/16 and a = 4, but 



= 0. 



• N on- Commutative Geometries: A = 2a ncg /E 2 and 
a — A, after renormalizing rn 9 and c. 

Of course, for different values of (A, a) we can parameter- 
ize other Lorentz-violating corrections to the dispersion 
relation. One might be naively tempted to think that 
a p 3, or p A correction to the above dispersion relation 
will induce a 1.5 or 2PN correction to the phase relative 
to the massive graviton term. This, however, would be 
clearly wrong, as p is the graviton's momentum, not the 
momentum of the members of a binary system. 

With this modified dispersion relation the modified 
graviton speed takes the form 
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AE a - 



To first order in A, this can be written as 



= 1 



2 4 

rrigC 



-AE° 



1 - 



m 2 c 4 



a/2 



(4) 



(5) 



In general, a modification of the dispersion relation will 
be accompanied by a change in either the Lorentz group 
or its action in real or momentum space. Lorentz- 
violating effects of this type are commonly found in quan- 
tum gravitational theories, including loop quantum grav- 
ity [29] and string theory [301 E] ■ 

Modifications to the standard dispersion relation are 
usually suppressed by the Planck scale, so one might 
wonder why one should study them. Recently, Collins, 
et al. [351 133] suggested that Lorentz violations in per- 
turbative quantum field theories could be dramatically 
enhanced when one regularizes and renormalizes them. 
This is because terms that would vanish upon renor- 
malization due to Lorentz invariance do not vanish in 
Lorentz-violating theories, leading to an enhancement af- 
ter renormalization |34j . 

Although this is an appealing argument, we prefer here 
to adopt a more agnostic viewpoint and simply ask the 
following question: What type of modifications would en- 
ter gravitational-wave observables because of a modified 
dispersion relation and to what extent can these devia- 
tions be observed or constrained by current and future 
gravitational-wave detectors? In view of this, we postu- 
late the parametrized dispersion relation of Eq. ([I]) . 

One can see that this model-independent dispersion 
relation can be easily mapped to all the ones described 
above, in the limit where E and p are large compared to 



and in the limit E m g it takes the form 



= 1 - 



E 2 



- AE C 



(6) 



Notice that if A > or if m 2 c 4 jE 2 > |A|£: Q ~ 2 , then the 
graviton travels slower than light speed. On the other 
hand, if A < and m 2 g c 4 /E 2 < |A|£ Q ~ 2 , then the gravi- 
ton would propagate faster than light speed. 



III. PROPAGATION OF GRAVITATIONAL 
WAVES WITH A MODIFIED DISPERSION 
RELATION 

We now consider the propagation of gravitational 
waves that satisfy the modified dispersion relation of 
Eq. ([T]). Consider the Friedman-Robertson- Walker back- 
ground 



ds 2 = -dt 2 + a 2 {t)[d X 2 + S 2 (x)(d6» 2 + i 



where a(t) is the scale factor with units of length, and 
£(x) is equal to %, sinx or sinhx if the universe is spa- 
tially flat, closed or open, respectively. Here and hence- 
forth, we use units with G = c = 1, where a useful con- 



version factor is 1M Q = 4.925 x 10 



1.4675 km. 
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In a cosmological background, we will assume that the 
modified dispersion relation takes the form 



A\p\ a 



(8) 



where \p\ = {gijP 1 ^) 1 ^ 2 ■ Consider a graviton emitted 
radially at \ = Xe an d received at x = 0. By virtue of 
the x independence of the t — x part of the metric, the 
component p x of its 4-momentum is constant along its 
worldline. Using E — p°, together with Eq. (jsj) and the 
relations 



p A 
~E 



dx 
dt ' 



p x = a 2 p x , 



we obtain 

dx = 
dt 



9 9 

m g a s 
Pi 



2-a 



(9) 



(10) 



where p 2 = a 2 (t e )(E 2 



A\p\"). The overall mi- 



nus sign in the above equation is included because the 
graviton travels from the source to the observer. 

Expanding to first order in (m g /E e ) <C 1, and 
A/p 2_a <C 1 and integrating from emission time (x = Xe) 
to arrival time (x = 0), we find 



dt 



1 



1 



te a(t) 2a\t e )El) t 

a-2 rt 



a(t)dt 



A a(t e )E c 



a(ty- a dt. (11) 



Consider gravitons emitted at two different times t e 
and t' e , with energies E e and E' e , and received at corre- 
sponding arrival times (Xe is t ne same for both) . Assum- 
ing At e = t e — t' e <C a/a, then 



At a = {l + Z) 



D„ 



At fi 



■Do 
2A 2 



1 

71 



i 

f 2 



2\i 



1; 



(12) 



where Z = aa/a(t e ) — 1 is the cosmological redshift, and 
where we have defined 



(13) 



and where m g /E e = (A g / e ) _1 , with f e the emitted 
gravitational- wave frequency, E e = hf e and A 9 = h/m g 
the graviton Compton wavelength. Notice that when 
a = 2, then the A correction vanishes. Notice also that 
Aa always has units of length, irrespective of the value of 
a. The distance measure D a is defined by 



D a = 



1 



rt(, 



l-Q 



a(t)^ a dt 



(14) 



where ao = a(t a ) is the present value of the scale factor. 
For a dark energy-matter dominated universe D a and the 
luminosity distance have the form 

Da - , d5) 



#o Jo ^/n M (i + z') 3 + n A 



D, 



l + Z 



dz' 



Ho Jo ^a/(1 + ^') 3 + ^a' 



(16) 



where Ho « 72 km s _1 Mpc - is the value of the Hubble 
parameter today and VLm = 0.3 and 51a = 0.7 are the 
matter and dark energy density parameters, respectively. 

Before proceeding, let us comment on the time shift 
found above in Eq. (12 1. First, notice that this equa- 
tion agrees with the results of [TT] in the limit A — > 0. 
Moreover, in the limit a ~ > 0, our results map to those 
of [TT] with the relation A^ 2 — > A^ 2 + A^ 2 . Second, no- 
tice that in the limit a — > 2, the (a(t e )E e ) 2 ~ a in Eq. ( 11 ) 
goes to unity and the A correction becomes frequency 
independent. This makes sense, since in that case the 
Lorentz-violating correction we have introduced acts as 
a renormalization factor for the speed of light. 



IV. MODIFIED WAVEFORM IN THE 
STATIONARY PHASE APPROXIMATION 

We consider the gravitational-wave signal generated 
by a non-spinning, quasi-circular inspiral in the post- 
Newtonian approximation. In this scheme, one assumes 
that orbital velocities are small compared to the speed of 
light (v <C 1) and gravity is weak (m/r <C 1). Neglecting 
any amplitude corrections (in the so-called restricted PN 
approximation) , the plus- and cross-polarizations of the 
metric perturbation can be represented as 



h(t) ee A(t)e- mt \ 

$(i) EE * C + 27T^ f(t)dt, 



(17) 
(18) 



where A(t) is an amplitude that depends on the 
gravitational-wave polarization (see e.g. Eq. (3.2) 
in [UJ), while f(t) is the observed gravitational- wave 
frequency, and $ c and t c are a fiducial phase and fidu- 
cial time, respectively, sometimes called the coalescence 
phase and time. 



The Fourier transform of Eq. (17 1 can be obtained an- 
alytically in the stationary-phase approximation, where 
we assume that the phase is changing much more rapidly 
than the amplitude [23 [36] • We then find 




»*(/) 



(19) 



where / is the gravitational-wave frequency at the detec- 
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tor and 

A(t) 



4 M c 



5 a S(K e ) 
rf 



*(/) = 2tt / (t - i c )4f + 2vr/t c 



(20) 
(21) 



In these equations, M. e = jj 3 ' 5 m is the c/wrp mass of the 
source, where 77 = m\mil{m\ + TO2) is the symmetric 
mass ratio. 



We can now substitute Eq. ( 12 ) into Eq. ( 21 ) to relate 



the time at the detector to that at the emitter. Assuming 
that q / 1, we find 

*aMf) = 2lT J f (te-t ec )df ( - 

1 ttLL 



(22) 



while for a 



(1 - a) f}- a \ 2 - a 
1, we find 



2nft c — $ c - 



7T 

4' 



*a=l(/) = 2^ / (t e -t ec )df, 



A A 



111 



A 



27T/i c - i c 



(23) 



7T 

4 



The quantities (i c ,i c ) and (<p c ,<p c ) are new coalescence 
times and phases, into which constants of integration 
have been absorbed. 

We can relate t e —t ec to f e by integrating the frequency 
chirp equation for non-spinning, quasi-circular inspirals 
from general relativity |llj : 



dfe 96 
dt e 5irM 2 



(*M e f B ) U ' 3 



1 



743 
336 



11 



+ 47r(7rm/ e 



r?)(7rm/ e ) 2 / 3 



(24) 



where we have kept terms up to 1PN order. In the cal- 
culations that follow, we actually account for corrections 
up to 3.5PN order, although we don't show these higher- 
order terms here (they can be found e.g. in [37]). 

After absorbing further constants of integration into 
(i c , <& c , i c , $ c ), dropping the bars, and re-expressing ev- 
erything in terms of the measured frequency / at the 
detector [note that f 1 ' 2 = (d/ e /rft e ) 1/2 /(l + Z)\, we ob- 
tain 

_/ i(/)e**(/\ for0</</„ 
0, for / > /„ 



Hf) 

with the definitions 

A{f) = eAu- v \ A 
*(/) = *qr(/) + **(/), 

*<»(/) - ^fU - $ 



/ max 
/ max i 



(25) 



7T 

4 



7T M 2 

-5/3 



(26) 



OC 

x i c « 

n=0 



128 

41n(u)]« n/3 > 



(27) 



where the numerical coefficient e = 1 for LIGO and 
ET, but e = y/3/2 for a LIS A- like mission (because 
when one angle-averages, the resulting geometric fac- 
tors depend slightly on the geometry of the detector). 
The coefficients (c n ,t n ) can be read up to n = 7 for 
example from Eq. (3.18) in [37]. In these equations, 
u = nA4f is a dimensionless frequency, while Ai is the 
measured chirp mass, related to the source chirp mass 
by M = (1 + Z)M e - The frequency f max represents 
an upper cut-off frequency where the PN approximation 
fails. 

The dephasing caused by the propagation effects takes 
a slightly different form depending on whether a =/= 1 or 
a = 1. In the general a ^ 1 case, we find 



where the parameters /3 and £ are given by 
ir 2 D M 

TT 2 - a D„ M 1 - 01 







/I 



(1 - a) \ 2 ~ a (1 + Z) 



\l-a 



In the special a = 1 case, we find 

<5*a=l(/) = -/3w _1 + Ca=lMu) , 

where /3 remains the same, while 

7r£>i 



A,: 



(28) 



(29) 



(30) 



(31) 



(32) 



and we have re-absorbed a factor into the phase of coa- 
lescence. 

As before, notice that in the limit A 



0, Eq. (28) 



reduces to the results of [TT] for a massive graviton. Also 
note that, as before, in the limit a — > 0, we can map our 
results to those of [11] with A^ 2 — > A^ 2 + A^ 2 , i.e. in 
this limit, the mass of the graviton and the Lorentz- 
violating A term become 100% degenerate. In the limit 
a — > 2, Eq. (12 1 becomes frequency-independent, which 
then implies that its integral, Eq. (21), becomes linear 
in frequency, which is consistent with the a — > 2 limit 
of Eq. ( 28 1 . Such a linear term in the gravitational- wave 



phase can be reabsorbed through a redefinition of the 
time of coalescence, and thus is not observable. This is 
consistent with the observation that the dispersion rela- 
tion with a = 2 is equivalent to the standard massive 
graviton one with a renormalization of the speed of light. 
When a = 1, Eq. |l2"l ) leads to a 1/f term, whose integral 
in Eq. ( 21 ) leads to a ln(/) term, as shown in Eq. ( 23 ) . Fi- 



nally, notice that, in comparision with the phasing terms 
that arise in the PN approximation to standard general 
relativity, these corrections are effectively of (l+3a/2)PN 
order, which implies that the a = term leads to a 1PN 
correction as in [TT], the a — 1 case leads to a 2.5PN 
correction, the a — 3 case leads to a 5.5PN correction 
and a = 4 leads to a 7PN correction. This suggests that 
the accuracy to constrain Aa will deteriorate very rapidly 
as a increases. 
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Connection to the Post-Einsteinian Framework 

Recently, there has been an effort to develop a 
framework suitable for testing for deviations from gen- 
eral relativity in gravitational-wave data. In analogy 
with the parametrized post-Newtonian (ppN) frame- 
work [TU1 [3HH33 > the parametrized post-Einsteinian 
(ppE) framework |12l [13l |43] suggests that we deform 
the gravitational-wave observable away from our GR ex- 
pectations in a well-motivated, parametrized fashion. In 
terms of the Fourier transform of the waveform observ- 
able in the SPA, the simplest ppE meta-waveform is 

KM = (1 + a ppE U^) e <*GR(/)+iA,pE « 6 - E f 

(33) 

where (a ppE , a ppE , /3 ppE , 6 ppE ) are ppE, theory parameters. 
Notice that in the limit a ppE — > or /3 ppE — > 0, the ppE 
waveform reduces exactly to the SPA GR waveform. The 
proposal is then to match-filter with template families of 
this type and allow the data to select the best-fit ppE pa- 
rameters to determine whether they are consistent with 
GR. 

We can now map the ppE parameters to those obtained 
from a generalized, Lorentz- violating dispersion relation: 

a ppE = ,3 ppE = -C & ppE = a — 1 . (34) 

Quantum-gravity inspired Lorentz- violating theories sug- 
gest modified dispersion exponents a — 3 or 4, to lead- 
ing order in E/m g , which then implies ppE parameters 
6 ppE = 2 and 3. Therefore, if after a gravitational wave 
has been detected, a Bayesian analysis with ppE tem- 
plates is performed that leads to values of 6 ppE that peak 
around 2 or 3, this would indicate the possible presence 
of Lorentz violation [13] . Notice however that the a = 1 
case cannot be recovered by the ppE formalism without 
generalizing it to include lnw terms. Such effects are 
analogous to memory corrections in PN theory. 

At this point, we must spell out an important caveat. 
The values of a that represent Lorentz violation for 
quantum-inspired theories (a = 3,4) correspond to very 
high PN order effects, i.e. a relative 5.5 or 7 PN correc- 
tion respectively. Any gravitational-wave test of Lorentz 
violation that wishes to constrain such steep momen- 
tum dependence would require a very accurate (high PN 
order) modeling of the general relativistic waveform it- 
self. In the next section, we will employ 3.5 PN accu- 
rate waveforms, which are the highest-order known, and 
then ask how well £ and /3 can be constrained. Since 
we are neglecting higher than 3.5 PN order terms in the 
template waveforms, we are neglecting also any possi- 
ble correlations or degeneracies between these terms and 
the Lorentz-violating terms. Therefore, any estimates 
made in the next section are at best optimistic bounds 
on how well gravitational-wave measurements could con- 
strain Lorentz violation. 



V. CONSTRAINING A MODIFIED GRAVITON 
DISPERSION RELATION 

In this section, we perform a simplified Fisher analysis, 
following the method outlined for compact binary inspi- 
ral in |4*4Tl4*^] , to get a sense of the bounds one could place 
on (A 9 ,Aa) given a gravitational- wave detection that is 
consistent with general relativity. We begin by summa- 
rizing some of the basic ideas behind a Fisher analysis, 
introducing some notation. We then apply this analysis 
to an Adv. LIGO detector, an ET detector and a LISA- 
like mission. 



A. General Considerations 

Given a noise power spectrum, S n (f), we can define 
the inner product of signals hi and h 2 &s 

{h 1 \h 2 ) = 2j^ — df, (35) 

where hi and hi are the Fourier transforms of signals 1 
and 2 respectively and star superscript stands for com- 
plex conjugation. The SNR for a given signal h is simply 

p[h] = (h\h) l /\ (36) 

If the signal depends on a set of parameters 9 a that we 
wish to estimate via matched filtering, then the root- 
mean-square error on parameter 9 a in the limit of large 
SNR is (no summation over a implied here) 

A9 a = y/({8 a - (0 a )) 2 ) = . (37) 

The quantity S aa is the (a, a) component of the variance- 
covariance matrix, which is the inverse of the Fisher in- 
formation matrix, r a ;,, defined as 



/ dh 


1 dh \ 







(38) 



The off-diagonal elements of the variance-covariance ma- 
trix give the parameter correlation coefficients, which we 
define as 

We will work with an angle-averaged response func- 
tion, so that the templates depend only on the parame- 
ters: 

6 a = {XnA^ c ,htcMMM V ,(i,0 , (40) 

where each component of the vector 6 a is dimensionless. 
We recall that A is an overall amplitude that contains in- 
formation about the gravitational-wave polarization and 
the beam-pattern function angles. The quantities $ c and 
t c are the phase and time of coalescence, where /o is a fre- 
quency characteristic of the detector, typically a "knee" 



7 



frequency, or a frequency at which S n (f) is a minimum. 
The parameters M. and r\ are the chirp mass and sym- 
metric mass ratio, which characterize the compact binary 
system under consideration. The parameters (/3, £) de- 
scribe the massive graviton and Lorentz-violating terms 
respectively. 



The SNR for the templates in Eq. (251 is simply 

p = 2 e A (Mn)- 7/6 I^) 1 / 2 S^ 1 ^ 2 , (41) 

where we have defined the integrals 



-9/3 



9{x) 



(42) 



with x = f/fo- The quantity g(x) is the rescaled power 
spectral density, defined via g(x) = Sh(f)/So for the de- 
tector in question, and Sq is an overall constant. When 
computing the Fisher matrix, we will replace the ampli- 
tude A in favor of the SNR, using Eq. (Ell). This will then 



lead to bounds on (/3, £) that depend on the SNR and on 
a rescaled version of the moments J(q) = I(q)/I(7). 

In the next subsections, we will carry out the in- 
tegrals in Eq. ( 42 ) , but we will approximate the lim- 
its of integration by certain x m - m and x max 15J. The 
maximum frequency will be chosen to be the smaller 
of a certain instrumental maximum threshold frequency 
and that associated with a gravitational wave emit- 
ted by a particle in an innermost-stable circular or- 
bit (ISCO) around a Schwarzschild black hole (BH): 
/max = 6~ 3 / 2 7r _1 ?7 3 / 5 A^ _1 . The maximum instrumen- 
tal frequency will be chosen to be (10 5 ,10 3 ,1) Hz for 
Ad. LIGO, ET and LISA-like, respectively. The mini- 
mum frequency will be chosen to be the larger of a cer- 
tain instrumental minimum threshold frequency and, in 
the case of a space mission, the frequency associated with 
a gravitational wave emitted by a test-particle one year 
prior to reaching the ISCO. The minimum instrumen- 
tal frequency will be chosen to be (10, 1, 10~ 5 ) Hz for 
Ad. LIGO, ET and a LISA-like mission, respectively. 

Once the Fisher matrix has been calculated, we will 
invert it using a Cholesky decomposition to find the 
variance-covariance matrix, the diagonal components of 
which give us a measure of the accuracy to which pa- 
rameters could be constrained. Let us then define the 
upper bound we could place on (/3, () as A/3 = A 1 / 2 /^ 
and AC = A 1 / 2 / p, where A and A are numbers. Com- 



bining these definitions with Eqs. (29) and (30), we find, 
for a/1, the bounds: 



Xg > 



Ar 2 < 



pD M 7T 

V (1 + Z) AV4 ' 

ll-alA 1 ^ M a - X 
tt 2 - q D a p (l + Z)*- 1 ' 



(43) 
(44) 



Notice that the direction of the bound on Aa itself de- 
pends on whether a > 2 or a < 2; but because A = 



(Aa//i)" 2 , all cases yield an upper bound on A. For the 
case a = 1 , we find 



(45) 



In the remaining subsections, we set /3 = and £ = in 
all partial derivatives when computing the Fisher matrix, 
since we derive the error in estimating /3 and C about the 
nominal or a priori general relativity values, (/3, £) = 
(0,0). 



B. Detector Spectral Noise Densities 



We model the Ad. LIGO spectral noise density via 



Shif) 
So 



( 10 16-4 (a; / o -7.9) 2 + 2 4 x 1Q -62 x -50 

+0.08ar 4 - 69 



ioq qc| 1-0. 23x^+0 . 0764a: 4 
~' ) '' ' 1+0. 17a; 2 

CO, 



/ > f s ; 

I < fs, 



(46) 



Here, f = 215 Hz, S = 10~ 49 Hz" 1 , and f s = 10 Hz is 
a low- frequency cutoff below which Sh(f) can be consid- 
ered infinite for all practical purposes 

The initial ET design postulated the spectral noise 
density |4"T] 




where f = 100 Hz, S = 10- 50 Hz _1 , f s = 1 Hz, and 



ai = 


2.39 x 10~ 27 , 


bi 


= -15.64, 




0.2 


= 0.349, 


b 2 


= -2.145, 






= 1.76, 


b 3 


= -0.12, 




«4 


= 0.409, 


bi 


= 1.10. 


(48) 



min|^ NSA (/)/exp (-^ ssion dN/df) , 
5r A (/) + 5f(/)}+5r Sal (/)- (49) 



The classic LISA design had an approximate spectral 
noise density curve that could be modeled via (see eg. [13 

my 



s h (f) = 

where 

sr A (f) 



9.18 x 10~ 52 



+9.18 x 10 



-38 



1 Hz 

/ 
1 Hz 



1.59 x 10" 



Hz" 



Sf'U) = 2.1 x ID- 



S'' 



cx— gal 



(/) = 4.2 x 10 



-47 



/ 

1 Hz 

/ 
1 Hz 



-7/3 



Hz" 1 



-7/3 



Hz" 1 



(50) 
(51) 
(52) 




FIG. 1: ET (left) and LISA (right) spectral noise density curves for the classic design (dotted) and the new NGO design 
(solid). 



and 

with A/ = T~ ission the bin size of the discretely Fourier 
transformed data for a classic LISA mission lasting a time 
^mission and k ~ 4.5 the average number of frequency bins 
that are lost when each galactic binary is fitted out. 

Recently, the designs of LISA and ET have changed 
somewhat. The new spectral noise density curves can be 
computed numerically [4"9TIS"T] and are plotted in Fig. [I] 
Notice that the bucket of the NGO noise curve has shifted 
to higher frequency, while the new ET noise curve is more 
optimistic than the classic one at lower frequencies. The 
spikes in the latter are due to physical resonances, but 
these will not affect the analysis. In the remainder of this 
paper, we will use the new ET and NGO noise curves to 
estimate parameters. 



C. Results 

We plot the bounds that can be placed on £ in Fig. [2] 
as a function of the a parameter. The left panel corre- 
sponds to the bounds placed with Ad. LIGO and p = 10 
(Dl ~ 160 Mpc, Z ~ 0.036 for a double neutron-star in- 
spiral), the middle panel corresponds to ET and p — 50 
(D L ~ 2000 Mpc, Z ~ 0.39 for a double 1OM BH inspi- 
ral) and the right panel corresponds to NGO and p = 100 
(D L ~ 20, 000 Mpc, Z ~ 2.5 for a double 10 5 M o BH in- 
spiral). When a = or a = 2, £ cannot be measured 
at all, as it becomes 100% correlated with either stan- 
dard massive graviton parameters. Thus we have drawn 
vertical lines in those cases. As the figure clearly shows, 
the accuracy to which £ can be measured deteriorates 



rapidly as a becomes larger. In fact, once a > 4, we 
find that £ cannot be confidently constrained anymore 
because the Fisher matrix becomes non-invertible (it's 
condition number exceeds 10 16 ). 

Attempting to constrain values of a > 5/3 becomes 
problematic not just from a data analysis point of view, 
but also from a fundamental one. The PN templates 
that we have constructed contain general relativity phase 
terms up to 3.5 PN order. Such terms scale as it 2 / 3 , which 
corresponds to a — 5/3. Therefore, trying to measure 
values of a > 5/3 without including the corresponding 
4PN and higher-PN order terms is not well-justified. We 
have done so here, neglecting any correlations between 
these higher order PN terms and the Lorentz-violating 
terms, in order to get a rough sense of how well Lorentz- 
violating modifications could be constrained. 

The bounds on f3 and £ are converted into a lower 
bound on X g and and upper bound on Aa in Table |TT] 
for a = 3 and binary systems with different component 
masses. Given a gravitational-wave detection consistent 
with general relativity, this table says that X g and Aa 
would have to be larger and smaller than the numbers in 
the seventh and eight columns of the table respectively. 
In addition, this table also shows the accuracy to which 
standard binary parameters could be measured, such as 
the time of coalescence, the chirp mass and the symmetric 
mass ratio, as well as the correlation coefficients between 
parameters. Different clusters of numbers correspond to 
constraints with Ad. LIGO (top), New ET (middle) and 
NGO (bottom - see caption for further details). 

Although Fig. [2] suggests bounds on ( of O(10 3 - 10 5 ) 
for the a = 3 case, the dimensional bounds in Table [TT] 
suggest a strong constraint on Aa- This is because in con- 
verting from £ to Aa one must divide by the D3 distance 
measure. This distance is comparable to (but smaller 
than) the luminosity distance, and thus, the longer the 
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10 M o +10 M D — i — 
10 M +100 M Q • 
- 100 M o +100 M Q * 














































New ET 



10* M. 
10^ M. 

10 6 M, 



1+10° M„ 
,+10^ M B 
,+10 b M n 



FIG. 2: Bounds on parameter £ for different values of a, using AdLIGO and p — 10 (left panel), ET and p = 50 (center panel), 
and NGO and p = 100 (right panel). Vertical lines at a = (0, 2) show where the £ correction becomes 100% degenerate with 
other parameters. Each panel contains several curves that show the bound for systems with different masses. 



Detector 


mi 


1712 


A(j> c 


At c 


AM/M 


A-q/n 


AA g 


AA A 


CMri 


cm 


C V 


C.MC 


c„ c 


C/3C 


Ad. LIGO 


1.4 


1.4 


3.61 


1.80 


0.0374% 


6.80% 


3.34 


0.911 


-0.962 


-0.991 


0.989 


-0.685 


0.803 


0.740 




1.4 


10 


3.34 


9.99 


0.267% 


12.8% 


3.48 


4.36 


-0.977 


-0.993 


0.917 


-0.830 


0.923 


0.875 




10 


10 


4.16 


31.0 


2.40% 


72.2% 


3.53 


8.40 


-0.978 


-0.994 


0.995 


-0.874 


0.947 


0.915 


ET 


10 


10 


0.528 


1.59 


0.0174% 


1.70% 


4.15 


0.0286 


-0.952 


-0.986 


0.988 


-0.742 


0.875 


0.813 




10 


100 


1.12 


44.5 


0.259% 


6.67% 


2.58 


1.38 


-0.974 


-0.993 


0.993 


-0.872 


0.951 


0.915 




100 


100 


5.23 


203 


4.03% 


67.6% 


1.86 


4.12 


-0.983 


-0.995 


0.996 


-0.914 


0.969 


0.947 


NGO 


10 4 


10 4 


0.264 


1.05 


0.00124% 


0.368% 


4.06 


0.266 


-0.957 


-0.990 


0.986 


-0.636 


0.761 


0.687 




10 4 


10 5 


0.264 


5.42 


0.00434% 


0.383% 


5.04 


1.81 


-0.955 


-0.991 


0.984 


-0.757 


0.884 


0.809 




10 5 


10 5 


0.295 


9.54 


0.0163% 


1.33% 


6.12 


2.33 


-0.944 


-0.983 


0.986 


-0.749 


0.891 


0.823 




10 5 


10 6 


0.351 


142 


0.0574% 


2.03% 


5.93 


118 


-0.961 


-0.990 


0.989 


-0.938 


0.942 


0.891 




10 6 


10 6 


0.415 


228 


0.138% 


5.33% 


8.30 


286 


-0.956 


-0.986 


0.990 


-0.820 


0.935 


0.885 



TABLE II: Root-mean-squared errors for source parameters, the corresponding bounds on A s and Aa, and the correlation 
coefficients, for the case a = 3 and for systems with different masses in units of Mq. The top cluster uses the Ad. LIGO S„(f), 
p = 10, A 9 is in units of 10 12 km, Aa is in units of 10~ 16 km and At c is in msecs. The middle cluster uses the ET S n (f), p = 50, 
X g is in units of 10 13 km, Aa is in units of 10 -15 km and At c is in msecs. The bottom cluster uses a NGO S„(f), p = 100, A 9 
is in units of 10 15 km, Aa is in units of 10 _10 km and At c is in sees. 



graviton propagates the more sensitive the constraints 
are to possible Lorentz violations. Second, notice that 
the accuracy to which many parameters can be deter- 
mined, e.g. (t c , AM, A77), degrades with total mass be- 
cause the number of observed gravitational-wave cycles 
decreases. Third, notice that the bound on the graviton 
Compton wavelength is not greatly affected by the inclu- 
sion of an additional parameter in the a = 3 case, and is 
comparable to the one obtained in [TT] for LIGO. In fact, 
we have checked that in the absence of Aa we recover 
Table II in [11]. 

We now consider how these bounds behave as a func- 
tion of the mass ratio. Figure [3] plots the bound on 
the graviton Compton wavelength (left) and the Lorentz- 
violating Compton wavelength Aa (right) as a function of 
r\ for Ad. LIGO and a — 3, with systems of different total 
mass. Notice that, in general, both bounds improve for 
comparable mass systems, even though the SNR is kept 
fixed. 

With all of this information at hand, it seems likely 



that gravitational-wave detection would provide useful 
information about Lorentz-violating graviton propaga- 
tion. For example, if a Bayesian analysis were carried out, 
once a gravitational wave is detected, and the ppE pa- 
rameters peaked around & ppE — 2 or 3, this could possibly 
indicate the presence of some degree of Lorentz violation. 
Complcmentarily, if no deviation from general relativity 
is observed, then one could constrain the magnitude of 
A to interesting levels, considering that no bounds exist 
to date. 



VI. CONCLUSIONS AND DISCUSSION 

We studied whether Lorentz symmetry-breaking in the 
propagation of gravitational waves could be measured 
with gravitational waves from non-spinning, compact bi- 
nary inspirals. We considered modifications to a massive 
graviton dispersion relation that scale as Ap a , where p is 
the graviton's momentum while (A, a) are phenomeno- 
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FIG. 3: Bounds on A 9 (left) and Aa (right) as a function of r) for different total masses, Ad. LIGO, p = 10 and a — 3. 



logical parameters. We found that such a modification 
introduces new terms in the gravitational-wave phase 
due to a delay in the propagation: waves emitted at 
low frequency, early in the inspiral, travel slightly slower 
than those emitted at high frequency later. This results 
in an offset in the relative arrival times at a detector, 
and thus, a frequency-dependent phase correction. We 
mapped these new gravitational-wave phase terms to the 
recently proposed ppE scheme, with ppE phase parame- 
ters 6 ppE = a — 1. 

We then carried out a simple Fisher analysis to 
get a sense of the accuracy to which such dispersion 
relation deviations could be measured with different 
gravitational-wave detectors. We found that indeed, 
both the mass of the graviton and additional dispersion 
relation deviations could be constrained. For values of 
a > 4, there is not enough information in the waveform 
to produce an invertible Fisher matrix. Certain values of 
a, like a — and 2, also cannot be measured, as they 
become 100% correlated with other system parameters. 

In deriving these bounds, we have made several ap- 
proximations that force us to consider them only as 
rough indicators that gravitational waves can be used 
to constrain generic Lorentz-violation in gravitational- 
wave propagation. For example, we have not accounted 
for precession or eccentricity in the orbits, the merger 
phase of the inspiral, the spins of the compact objects 
or carried out a Bayesian analysis. We expect the inclu- 
sion of these effects to modify and possibly worsen the 
bounds presented above by roughly an order of magni- 
tude, based on previous results for bounds on the mass of 
the graviton [TT1 UM [TBI I52"H55] . However, the detection 
of N gravitational waves would lead to a ^/~N improve- 
ment in the bounds 19J, while the modeling of only the 
Lorentz-violating term, without including the mass of the 
graviton, would also increase the accuracy to which Aa 



could me measured [13]. 

Future work could concentrate on carrying out a more 
detailed data analysis study, using Bayesian techniques. 
In particular, it would be interesting to compute the 
evidence for a general relativity model and a modified 
dispersion relation model, given a signal consistent with 
general relativity, to see the betting-odds of the signal 
favoring GR over the non-GR model. A similar study 
was already carried out in |13j . but there a single ppE 
parameter was considered. Another interesting avenue 
for future research would be to consider whether there 
are any theories (quantum-inspired or not) that predict 
fractional a powers or values of a different from 3 or 4. 
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